On the D-wave state component of the deuteron in 
the Nambu-Jona-Lasinio model of light nuclei 



A. N. Ivanov , V. A. Ivanova *, H. Oberhummer 
N. I. Troitskaya *, M. Faber ^ 

February 8, 2008 



Atominstitut der Osterreichischen Universitdten, Arbeitsbereich Kernphysik und 
Nukleare Astrophysik, Technische Universitdt Wien, 
Wiedner Hauptstr. 8-10, A-1040 Wien, Osterreich 

Abstract 
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is calculated in the Nambu-Jona-Lasinio model of light nuclei. For the ratio of D- 
to S-state deuteron wave functions we obtain equal to % = 0.0238. This agrees well 
with the phenomenological value % = 0.0256 ±0.0004 quoted by Kamionkowski and 
Bahcall (ApJ. 420, 884 (1994)). 
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1 Introduction 



The Nambu-Jona-Lasinio model of light nuclei or differently the nuclear Nambu- 
Jona-Lasinio (NNJL) model suggested in [1-3] represents a relativistically covariant quan- 
tum field theoretic approach to the description of low-energy properties and interactions 
of the deuteron and light nuclei. The NNJL model is fully motivated by QCD [1]. The 
deuteron appears in the nuclear phase of QCD as a neutron-proton collective excitation, 
the Cooper np-pair, induced by a phenomenological local four-nucleon interaction. The 
NNJL model describes low-energy nuclear forces in terms of one-nucleon loop exchanges 
providing a minimal transfer of nucleon flavours from initial to final nuclear states and 
accounting for contributions of nucleon-loop anomalies which are completely determined 
by one-nucleon loop diagrams. The dominance of contributions of nucleon-loop anoma- 
lies to effective Lagrangians of low-energy nuclear interactions is justified in the large Nc 
expansion, where Nq is the number of quark colours. 

Nowadays there is a consensus concerning the existence of non-nucleonic degrees of 
freedom in nuclei [4] . The non-nucleonic degrees of freedom can be described either within 
QCD in terms of quarks and gluons [5] or in terms of mesons and nucleon resonances [6]. 
In the NNJL model the non-nucleonic degrees of freedom of nuclei have been investigated 
in terms of the A(1232) resonance and calculated the contribution of the AA component 
to the deuteron [2]. The obtained result P(AA) = 0.3 % agrees well with the experimental 
upper bound P(AA) < 0.4% [7] and other theoretical estimates [4]. 

As has been shown in [3] the NNJL model describes well low-energy nuclear forces for 
electromagnetic and weak nuclear reactions with the deuteron of astrophysical interest 
such as the neutron-proton radiative capture n + p — ► D + 7, the solar proton burning 
p + p^D+e + + z/ e , the pep-process p + e~+p^D + z/ e and reactions of the 
disintegration of the deuteron by neutrinos and anti-neutrinos caused by charged u e + D 
— >e~ + p + p, z/ e + D^e + + n-|-n and neutral z/ e (z/ e ) + D — > u e (P e ) + n + p weak 
currents. 

The important problem which has not been jet clarified in the NNJL model is related 
to the value of the contribution of the D-wave state to the wave function of the deuteron. 
In this paper we fill this blank. In section 2 we calculate the contribution of the D- 
wave state to the wave function of the deuteron. We use a relativistically covariant 
partial-wave analysis developed by Anisovich et al. [8] for the description of nucleon- 
nucleon scattering. The fraction of the D-wave state of the deuteron wave function 
relative to the S-wave one we obtain equal to % = 0.0238. This agrees well with the value 
rid = 0.0256 ± 0.0004 quoted by Kamionkowski and Bahcall [9]Q who used this parameter 
for the phenomenological description of the realistic wave function of the deuteron in 
connection with the calculation of the astrophysical factor 5 , pp (0) for the solar proton 
burning p + p — ► D + e + + u e in the potential model approach. In the Conclusion we 
discuss the obtained result. 



1 The value rjd = 0.0256 ± 0.0004 was taken by Kamionkowski and Bahcall from Ref. [10]. 
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2 The D-wave state component of the deuteron 



The calculation of the value of the D-wave state contribution to the wave function of 
the deuteron we would carry out in terms of the amplitude of the transition n + p — > D. 
We show that the neutron-proton pair couples to the deuteron in both the S-wave state 
and the D-wave state with the fraction of the D-wave state agreeing with low-energy 
nuclear phenomenology. 

In the NNJL model the phenomenological Lagrangian of the npD interaction is defined 

by[i] 



--iipD 



[x) = —igy[p c (x)^n{x) — n c {x)^p{x)]D^(x) 



9t 



\p c (x)a^n(x) - n c {x)a^p{x)]Dl u {x) + h.c. (2.1) 



2 .Mm V ' V ' V ' ^ /J 



where D' (x), n(x) and p(x) are the interpolating fields of the deuteron, the neutron and 
the proton, {x) = d^D^x) — d v D^ (x) is the deuteron field strength. The phenomeno- 
logical coupling constant gy is related to the electric quadrupole moment of the deuteron 
Q D = 0.286 fm 2 , g\ = 2ir 2 Q D M^ [1], where M N = 940 MeV is the nucleon mass. The 
coupling constants gy and g T are connected by the relation [1] 

ST = ^f<?v, (2.2) 

which is valid at leading order in the large Nc expansion [1]. 
The amplitude of the transition n + p — > D is determined by 

Ih \ \r <(\\\h rr-h rr\ M M k n, ^n) + PQ%, <Tp) ~> D(fc D , Ap)) 

{Kb, A d \L np i){v) \k p , a p , fc n , cr n ; = , , (2.3) 

^2E D V2E n V2E p V 

where (.Ed, &d, Ad), (E p ,k p ,a p ) and (E n ,k n ,a n ) are energies, 4-momenta and polariza- 
tions of the deuteron, the proton and the neutron, respectively, V is a normalization 
space volume. The wave functions of the initial and the final states of the transition n + 
p — > D are given by 

\k p ,a p ; k n ,a n ) = a P (fcp, a p ) al(k n , <7 n )|0), 

(kn, Ad | = (0|a D (fc D ,A D ), (2.4) 

where a p (k p ,a p ) and a^(A; n , <r n ) are creation operators of the proton and the neutron, 
od(^d, Ad) is the annihilation operator of the deuteron and |0) is a vacuum wave function. 
The relativistically invariant amplitude M(n(k n ,<7 n ) +p(k p ,o- p ) — > D(&d, Ad)) reads 

M(n(k n , <r n ) + p(k p , a p ) -> D(/c D , A D )) = e*"(fc D , A D ) 

x ^2ig Y [u c (k n , a n )>y u u(k p , a p )\ - ^^["^n, o^(J^u(k p , a p )\ (k n + /c p ) M |, (2.5) 

where u c (k Q , <r n ) and w(A; p , cr p ) are bispinorial wave functions of the neutron and the proton 
with 4-momenta k n , k p and polarizations <j n , cr p ; e*^(/cD, Ad) is a 4-vector of polarization 
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of the deuteron with a 4-momentum kn and polarization Ad- The 4-momenta kn, k n and 
k v are related by kn = k n + k p due to conservation of energy and momentum. 

As has been shown by Anisovich et. al. [8] for neutron-proton scattering the neutron- 
proton densities describing the S- and D-wave states of a neutron-proton pair are equal 
to 



tf„( 3 Di;<r n ,CT p ) 



[u c (k n ,a n )S u u(k p ,a p )], 
[u c (k n ,a n )V u u(k p ,a p )}, 



(2.6) 



where S u and T> u are relativistically covariant operators of the projection onto the S-wave 
and the D-wave state, respectively [8]: 



2,k u 



2M N + y/s~ 



,3/2 



- (AM 2 - S )rf - (M N + y/s) h 



Here P = k p + k n , k = \ (k p — k a ), s = P 2 , P ■ k = and 



(2.7) 



X f) v 

lv =lu-P—- 

S 



(2.8) 



The neutron-proton densities Eq.( [2.6|) are normalized by the condition [8] 



i J ti{L^k p + M N )L»(-k n + M N )} (2vr)V 4) (P - k p - k» 



d 3 K 



d 3 k n 



(2tt) 3 2£ p (27r) 3 2E n 



Pl{s) 



(2.9) 



where = or X> M , the factor 3 in the denominator of the l.h.s. describes the number 
of the states of a neutron-proton density with a total momentum J = 1, 2 J + 1 = 3, and 
p,s(s) and pv{s) amount to 



87T 



Pv{s) 



1 (s-AM 2 \ 1/2 



s -AM, 



2\5/2 
N \ 



8vr \ s / 

In the center of mass frame of the neutron-proton pair the densities Eq. 
* ( 3 Si;<T n ,Op) 



(2.10) 
are equal to 



pi- 



= [u c (k n , a n )S u(k p , cr p )] = 0, 

^( 3 Si;cr n ,Op) = [u c (k Q ,a n )Su(k p ,a p )] = -j=<fl(a n )aip p (a 

$o( 3 Di;a n ,(Tp) = [u c (k n , a n )V u(k p , a p )) = 0, 

^( 3 D i;( T n ,(T p ) = [u c (k Q ,a n )Vu(k p ,a p )] = -^i(a n )[3 (a ■ v) v ~ v 2 a] ^ p (a p ) , (2.11) 



where v = k/y k 2 + = J 1 — 4M^/s and k are a relative velocity and a 3-momentum 
of the neutron-proton pair, <^ n (cr n ) and <£>p(cr p ) are spinorial wave functions of the neu- 
tron and the proton, respectively. It is obvious that the densities Eq.( [2.11|) describe the 
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neutron-proton pair in the S- and D-wave states with a total spin S = 1 and a total 
momentum J = 1. 

The neutron-proton densities Eq. fl2.llp are normalized by 

\ E E ^Sl^n^p)-^!^^) = 1, 

o- n =±l/2 Op=±l/2 

^ E E ^( 3 Di;a n ,a p ).^( 3 D i; a n ,a p ) = , 4 Jl-^j. (2.12) 

cr n =±l/2 cr p =±l/2 V S / 



The decomposition of the neutron-proton densities in the amplitude Eq.( |2.5| ) into the 
densities with a certain orbital momentum we would carry out at leading order in the 
large Nc expansion [1-3]. This would allow to consider the neutron and the proton as 
free particles obeying free equations of motion 

u c {k Q ,a n ){k n + M N ) = 0, 

(k p -M N )u(k p ,a p ) = 0. (2.13) 



In order to express the neutron-proton densities in the amplitude Eq.( |2.5| ) in terms of the 
projection operators Eq. (|2.7|) , first, we have to exclude the term containing a^. This can 
be carried out by using Gordon's identity 

(fa _|_ & )/■* k 
[u c (k n , a^cr^uikp, a p )} n p = -[u c (k n , a n )-i u u(k p) a p )} + — p- [u c (k n , a n )u(k p , a p )]. 

(2.14) 

Substituting Eq.(gTTJ) in Eq.(p75|) we get 

M(n(k n , a a ) + p(k p , (x p ) -> D(A; D , A D )) = 2% v + 2# T ) e^(/c D , A D ) 

x I [w c (A; n , cr n )j v u{k p , a p )} [u c (k Q , a n )u(k p , cr p )] i . (2.15) 

t SV + % J 

In terms of <Sj, and X>„ vectors 7^ and k u are determined by 

*" = ^ { . iM i-'>) s --\ sV - < 2 ' 16 ) 

Substituting Eq.( [2.16| ) in Eq. (|2.15|) and taking into account that [u c (k n , cr n )Pu(k p , a p )\ = 
we obtain 

M(n(k n , a n ) + p(k p , a p ) -> D(Ajd, A D )) = ^i{gy + 2^ T ) M N e*"(A; D , A D ) 
x {[w c (£; n , a n )S v u(k p , a p )] + ?? d [w c (£: n , a n )V„u(k p , a p )}} = 

= iV2i(gy + 2g T ) M N e* u (k D , A D ) [tt„( 3 Si; a n , a p ) + r/ d #„( 3 D i; a n , a p )}, (2.17) 

where rjd describes the fraction of the D-wave state in the wave function of the deuteron. 
It is equal to 

1 2g T -g Y 

Vd = — 7= 7; • (2.18) 
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For the derivation of Eqs.( [2.17 ) and (p. 18 ) we have set s = M§ and neglected the con- 
tribution of the binding energy of the deuteron in comparison with a nucleon mass M^. 
This means that 4M^- — s = when compared with M^. 

Using the relation Eq.(|2.2|) the parameter rja takes the value 



1 V3-V% 

571 7177!^ oom (219) 

This agrees well with the value rjd = 0.0256 ± 0.0004 that was used in low-energy nuclear 
phenomenology for the description of the realistic wave function of the deuteron within 
the potential model approach [9, 10]. 



3 Conclusion 

We have shown that the NNJL model describes well in agreement with low-energy 
nuclear phenomenology [10] such a fine structure of the deuteron as a contribution of the 
D-wave state. The calculation of the fraction of the D-wave state to the wave function of 
the deuteron we have carried out at leading order in the large Nq expansion [1] . This has 
allowed to treat the neutron and the proton as free particles on-mass shell [2] obeying 
free equations of motion. To the decomposition of the amplitude of the transition n 
+ p — ► D into the neutron-proton quantum field configurations having certain orbital 
momenta and corresponding to the S- and D-wave states, respectively, we have applied 
a relativistically covariant partial-wave analysis invented by Anisovich et al. [8] for the 
description of nucleon-nucleon scattering with nucleon-nucleon pairs coupled in the states 
with certain orbital momenta. 

The theoretical value of the D-wave state fraction in the wave function of the deuteron 
r]d = 0.0238 calculated in the NNJL model agrees well with low-energy nuclear phe- 
nomenology giving rjd = 0.0256 ± 0.0004 [10]. The former was quoted by Kamionkowski 
and Bahcall [9] for the parameterization of the realistic wave function of the deuteron 
in connection with the calculation of the astrophysical factor S pp (0) for the solar pro- 
ton burning p + p^D + e + + z/ e . The calculation of the contribution of the D-wave 
state fraction of the wave function of the deuteron in agreement with low-energy nuclear 
phenomenology testifies that the NNJL model describes to full extent low-energy tensor 
nuclear forces playing an important role in low-energy nuclear physics on the whole and 
for the existence of the deuteron, in particular [11]. 
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